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PART  I 


INTRODUCTION 


The  study  of  regional  wave  propagation  allows  to  know  the 
general  lineaments  of  the  internal  earth  structure:  remarkable 
discontinuities,  mantle  and  crust  dynamics  (Dziewonski,  1984; 
Jordan,  1981;  Jordan  et  al.  ,  1975,  1989;  Johnson,  1967,  1969; 
Burdich  et  al.,  1978;  Lay  et  al.,  1984). 

On  the  othter  hand  the  coda  of  P-waves  for  teleseisms 
(distance  larger  than  10®)  helps  to  establish  low  velocity  zones 
(LVZ)  beneath  different  structures:  those  LVZ  are  not  apparent,  or 
do  not  exist  at  all,  in  shield  structures  and  on  the  contrary  they 
are  ticker  in  tectonically  active  regions  (Helmberger,  197  3)  .  This 
data  informs  us  about  thermal  variations  related  to  lateral 
variations  in  the  upper  mantle. 

Regional  waves  (distance  of  10®  or  less)  help:  to  detect 
variations  in  the  crust;  more  precise  locations  of  seismic  events. 
r,ocal  phases  allow  to  know  with  more  detail  characteristics  close 
fo  the  station  or  to  the  station  network. 
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other  important  phases  independent  of  the  distance  (that  is, 
for  earthquakes,  either  local,  or  regional,  or  teleseism)  are 
identified  with  those  short  period  waves  Li,  Lg,  Rg,  what  are 
characteristics  of  continental  regions;  they  allow  a  better 
knowledge  of  low  velocity  layers  in  the  crust.  (Nuttli,  1986;  Payo 
Subiza,  1974;  Pomeroy  et  al.,  1980;  Alcocer,  1989;  Ayala,  1989; 
among  others) . 

Nevertheless,  those  phases  have  different  characteristics  for 
(-ifferent  complex  zones,  meanwhile  they  are  similar  for  more 
homogeneous  zones. 

South  America,  especially  the  Pacific  Coast,  is  really  complex 
when  considering  tectonic  structures,  crustal  thickness  and  Nazca 
plate  subduction.  Let  us  distinguish  two  types  of  geological 
structures:  stable  zones  (Guyana,  Brazil  and  Patagonia  Shields  and 
the  Altiplano) ;  strong  activity  zones  (Western  and  Eastern 
Cordil leras)  . 

Crustal  thickness;  40  km  at  the  north,  increases  to  65  km  in 
central  part  and  then  decreases  to  3  5  km  in  the  south  (Meissner  et 
al.,  1976;  Couch  et  al.,  1981;  James,  1971;  Wigger  et  al.,  1990; 
Fernandez  ,  1965;  Valez,  1982). 

Nazca  subducted  plate:  Subduction  dipping  almost  null  north 
;f  latitude  13°-  I5‘’s,  changes  to  35°  at  latitude  20°S;  dipping  is 
maximum  for  maximum  depth  between  500  to  700  km  (Deza,  1972; 
James,  1990;  Minaya,  1973). 
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Those  aspect  of  South  America  and  Pacific  Coast  complexity 
have  a  strong  influence  on  the  waves  recorded  in  LPB  station.  So 
the  identification  of  those  different  phases  is  one  of  the  highest 
importance . 


DATA 


1000  Seismic  events  which  occurred  in  South  America  and 
Pacific  Coast,  recorded  in  the  LPB  station  (La  Paz  -  Bolivia,  WVJSSN 
type)  were  analysed,  after  excluding  other  about  1000  events. 
Epicentral  data  were  read  from  the  International  Seismological 
Centre  (ISC)  Bulletins  from  1974  through  1988. 

Different  phase  analysis  was  achieved  on  the  3  components  of 
analogical  records. 

Distance  of  teleseisms  reach  36.5°;  regional  seisms  are  at 
distances  between  1°  and  10°.  Magnitude  mb  is  larger  than  4.  Depth 
may  go  down  to  650  km.  For  local  erthquakes  a  small  error  in  depth 
influences  records  too  much;  so,  to  avoid  that  difficulty,  events 
of  the  seismic-geographical  region  of  Flinn  and  Engdahl  3-120 
(Bolivia)  were  excluded  from  our  study. 

Seismic  records  were  digitized,  using  a  time  window  of  25  to 
35  sec,  being  sample  interval  0.05  sec. 
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Flinn  and  Engdahl  seismic  regions  studied  are:  7,  8,  9,  35. 


METHODOLOGY 

1.  Seismic  phases,  amplitude  (mm)  and  period  were  read  on 
analogical  short  period  records;  precision  depends  on  the 
abruptness  of  the  initial  phase.  Particle  motion  helped  to  correct 
time  readings  for  the  phase  after  the  first  arrival.  Epicentral 
distance  (km)  divided  by  travel  time  (second)  was  called  apparent 
velocity. 


2.  Phase  identification:  Trying  to  identify  t.he  regional 
phases,  frequently  P  was  quite  clear;  tor  distance  Larger  than  10 
often  pP  and  sP  were  detected;  other  regional  phases  derived  from 
?  could  not  be  identified;  but  certainly  other  phases  are  present, 
what  were  called  PI  and  P2 . 

A  phase  type  S  may  be  Sn,  but  was  not  labeled  so,  considering 
that  Moho  depth  is  not  known  reliably. 

We  have  to  remark  that  Lg  phase  is  clearly  recorded  for 
earthquakes  north  of  LPB  station,  independently  of  depth;  south  of 
LP3  Rg  amplitude  may  equal  that  of  Lg,  or  be  lesser  or  also 
completely  absent. 
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Particle  motion  helped  to  identify  phases  close  to  first 
arrival  (Figs.  1,  2,  3). 

3.  Analysis  of  crustal  tickness,  geology  and  tectonics  in  the 
region  and  Nazca  plate  subduction: 

Looking  for  the  changes  in  crustal  tickness,  let  us  see  a 
profile  from  north  to  south  (Fig.  4) : 

Along  the  parallel  B.S^N  Meissner  et  al.  (1976)  found  that 
oceanic  crust  is  20  km  thick  close  to  the  coast,  and  maximum 
tickness  of  40  km  in  the  continent. 

At  14 °S  he  found  from  west  to  east  20  km,  then  a  maximum  of 
60  km,  decreasing  to  30  km  in  the  Subandean. 

Forth  et  al.  (1990)  in  a  profile  along  the  parallel  2l'’s  find 
a  thickness  of  65  km  between  66*’  and  68.25°W;  at  the  coast  39  km; 
at  the  west  of  66 “w  30  km. 

Couch  et  al.  (1981)  calculate  a  thickness  of  45  km  in  the 
continent  at  latitude  26 °S  and  only  8  km  in  the  oceanic  crust. 

At  the  latitude  38 "S  the  same  author  calculates  a  maximum 
thickness  of  35  km  in  the  continent. 
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Geological  and  tectonic  changes  are  more  relevant  from  east 
to  west  than  from  north  to  south:  Shield  zones,  Subandean, 
Cordilleras  (Western,  Central  and  Eastern) ,  Coast  Cordillera,  coast 
and  transition  zone  from  continental  to  oceanic  crust. 

According  to  Deza  (1972),  Nazca  plate  has  three  transition 
zones:  between  0““  and  l‘’S;  between  13°  and  14 °S;  between  26°  and 
27 °S.  James  (1990)  considers  the  transition  to  begin  at  15 °S.  Those 
transition  zones  coincide  with  the  change  of  subduction  angle  of 
Nazca  plate. 

That  analysis  recommended  to  divide  seismic  regions  (Fig.  2) : 


SEISMIC  REGION 

ZONE 

GEOGRAPHIC  REGIONS 

7 

7a 

98,  partial  both  97  and  101 

7b 

96,  99,  100,  partial  both  97  and  101 

3 

3a 

’02,  104,  105,  108,  109 

3b 

103,  106,  107,  110,  111,  partial  both  112  and  116 

3c 

113,  119,  '20 

3d 

partial  '12,  114,  115,  partial  116,  117,  113 

8e 

122,  123,  124,  125,  partial  127,  123,  129 

3f 

partial  127,  130,  131,  132,  partial  136  through  14 

9 

9 

143  through  146 

55 

55 

528  through  531 

6 
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Earthquakes  located  in  seismic  zones  8a  and  35  have  no  phase 
derived  from  P.  Seismic  record  is  of  small  amplitude.  S  phase  is 
not  visible.  Lg  is  of  small  amplitude.  We  have  to  remar.k  that 
earthquakes  located  in  8a  have  a  path  to  LPB  different  of  those  in 
35.  Those  in  8a  are  originated  beneath  the  ocean  in  the  coast,  or 
in  the  subducted  plate;  those  in  zone  35  are  shield  surface  quakes. 

Other  phases  close  to  P  were  observed  for  the  earthquakes 
originated  in  all  the  other  provinces:  7a,  7b,  3b,  8c,  8d,  3e,  3f, 
3g  and  9.  But  there  are  so  many  differences  in  P  derived  phases 
that  they  claim  for  a  separate  description. 

Seismic  Region  7a.  After  P  a  larger  phase  arrives  with  an  interval 
corresponding  to  depth;  that  means  that  if  is  a  pP  (Fig.  6.1). 

See  the  gradation  of  amplitudes  P<pP<Lg  +  Rg  (both  later 
phases  indepent  from  depth) .  S<<P  .  For  the  corner  northeast  of 
this  region  P  is  fairly  large;  no  other  P  derived  phase  is 
observed;  Lg<P.  80%  Of  path  for  those  earthquakes  to  LPB  is  across 
shield  zones. 

Seismic  region  7b.  The  main  difference  with  to  appears  in  short 
period  surface  waves:  Surface  earthquakes  located  between  72°  and 
o9°W  and  7°  to  9°N  have  them  larger  than  pP  (Fig.  6.2). 
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Intermediate  earthquakes  have  smaller  short  period  surface 
waves . 

Seismic  region  8a.  Generally  no  P  derived  phase;  in  a  few  cases  a 
second  phase,  3  to  8  seconds  after  P  of  unknown  process  of 
generation. 

Lg  and  Rg  <<  P,  ;  S  very  small. 

Seismic  region  3b.  In  most  of  the  records  pP  is  observed;  in  some 
cases  another  phase  P,  between  P  and  pP,  5  to  9  seconds  after  P. 

Short  period  surface  waves,  especially  Rg,  are  clear  and 
large;  they  appear  for  intermediate  depth  as  well  as  for  surface 
quakes . 

Three  phases  after  P  are  observed,  but  in  several  cases  only 
the  phase  P  is  observed.  The  whole  record  has  a  small  amplitude 
(Fig.  6.3)  . 

Seismic  region  8c.  Most  of  the  earthquakes  of  this  region  are  deep. 
The  surface  or  intermediate  quakes  have  three  P  phases  (Fig.  6.6), 
being  the  second  one  pP  larger  than  the  other  two  (1st  and  3rd  this 
one  not  identified) .  S  phase  is  clear  for  deep  earthquakes  (more 
than  500  km),  but  emergent  for  surface  and  intermediate  quakes. 
Short  period  surface  waves  are  Larger  for  surface  and  intermediate 
quakes,  not  existent  for  deep  quakes. 
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Seismic  region  8d.  It  is  considered  especially  complex,  since 
record  amplitude  may  be  larger  or  small  for  the  same  magnitude;  ? 
derived  phases  may  be  clear  or  completelly  absent;  in  some  cases 
a  regular  pP  was  identified.  We  have  to  remark  that  Natca  plate 
dip  changes  gradually  within  this  region  from  horizontal  to  35'' 
(Fig.  6.5,8,17)  . 

Seismic  region  8e.  Records  show  two  phases  after  P,  one  P.  to  5 
seconds,  the  other  14  to  20  seconds  after  P.  Looking  to  the 
envelope  curve,  two  are  clearly  different: 

1  The  third  arrival  is  the  largest  one  (Fig.  6.12  through  16). 

2  The  first  P  arrival  is  the  largest  one  (Fig.  6.18). 

Short  period  surface  waves  are  smaller  than  S;  records  are  not 
clear,  especially  for  surface  waves,  without  doubt  because,  being 
regional  events,  wave  trains  recorded  analogically  may  not  be 
decompressed . 

Seismic  region  8f.  A  second  arrival  2  to  8  seconds  after  P  is 
independent  from  the  focal  depth.  Deep  earthquakes  (h>500  km;  do 
not  show  any  duplication  of  P,  that  is  to  say,  no  pP. 

Generally,  but  no  always,  the  first  phase  is  smaller  than  the 
second  one.  Short  period  surface  waves  Lg,  Rg  are  larger  than 
previous  phases  for  surface  quakes,  are  smaller  for  intermediate, 
are  absent  for  deep  ones.  No  Li  recorded  for  this  region.  (Fig.  6.9 
through  11) . 
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Seismic  region  3q.  Ocean  and  coast.  Generally  the  whole  record  is 
of  small  amplitude.  A  second  phase,  larger  than  P,  corresponds  to 
pP  (Fig.  6.8). 

Seismic  region  9.  3  To  9  seconds  after  P,  another  phase  is 
observed,  and  a  third  one  after  18  seconds.  Generally  the  coda  is 
short,  signal  amplitude  small.  Short  period  surface  waves  are  not 
observed  at  all. 

Seismic  region  35.  This  geological  zone  is  quiet,  well  known, 
consisting  mostly  in  Guyana  and  Brazil  shields.  One  P  phase  of 
small  amplitude  is  observed  (Fig.  6.4). 

Phases  Li,  Lg  and  Rg  are  clear.  S  phase  has  small  amplitude 
and  a  rather  long  period,  so  it  is  read  on  long  period  records. 
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CONCLUSIONS  AND  INTERPRETATION 

Analysis  of  the  phases  following  the  first  P  arrival  for 
earthquakes  occurring  beneath  South  America  and  Pacific  coast 
recorded  at  LPB  station  has  shown  that  these  phases  are  related  to 
the  South  American  and  Nazca  plate  complex  structure. 

The  complexity  of  Andean  Cordillera  explains  the  strong 
attenuation  of  waves  transmitted  along  a  path  parallel  to  its  axis 
across  transition  zones.  Variation  of  crustal  thickness  appears  to 
influence  more  seismic  waves  than  source  or  station  conditions. 

Problably  P,  has  its  origin  in  a  canalization  of  waves  in  the 
subducted  Nazca  plate;  also  a  low  velocity  layer  in  the  upper 
mantle  was  suggested  for  teleseisms  of  intermediate  dept.h.  For 
regional  earthquakes  of  surface  focus  rather  reflexion  in  the 
multiply  layered  continental  crust  would  be  the  responsible  of 
those  phases;  that  would  also  explain  the  delay  in  many  earthquakes 
recorded  with  small  amplitude. 

P,  possibly  is  originated  in  the  S  to  P  conversion;  but  this 
possibility  does  not  explain  P^  presence  for  deep  earthquakes. 
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In  a  few  earthquakes  P  is  faster  than  expected  according  the 
tables.  These  cases  have  an  irregular  geographical  distribution  net 
coincident  with  the  division  of  South  America  related  to  other 
seismic  aspects.  Possibly  a  high  velocity  layer  in  the  mantle  may 
explain  that  (20°  discontinuity  400  km  deep.  Lehmann,  1970;  Julian 
et  al.,  1968;  Fahmi,  1964). 

Nazea  plate  geometry  is  important  related  to  conversion  of 
wave  type  or  its  reflexion,  mainly  if  focus  is  at  an  intermediate 
depth. 

On  the  other  hand  the  wedge  of  the  mantle  between  the  Moho 
and  the  Benioff  zone  may  disturb  seismic  ray  propagation  for  deeper 
events,  those  originated  in  plate  contact. 

The  small  amplitude  and  the  absence  of  other  P  derived  phases 
in  the  region  a  possibly  are  explained  by  a  strong  absorption  of 
energy  in  the  transition  zone  oceanic  to  continental  crust. 

The  region  35  produces  very  small  P,  without  any  acceptable 
explanation  for  the  moment. 

See  in  a  tabular  form  the  behavior  of  South  American  regions 
concerning  P  derived  phases  at  LPB;  those  of  unknown  origin  will 
be  called  P,  and  P, . 
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ZONE 

One  P  arrival  8a 

35 

Two  P  phases  7a 

7b 

Three  P  phases  8b 

8c 

8d 

8e 

8f 

8g 
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P 

P 

P 

P 

P 

P 

P 

P 

P 

P 
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P-PHASES 


Pi 


Pi 

Pi 

Pl 


Pi 


pP 

pP 

pP 

PP 

pP 

pP 

pP 


P^ 
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A  paper  was  presented  in  the  General  Assembly  XX  of  the  lUGG 
in  3  poster  session  of  August  22,  whose  abstract  is  as  follows: 

REGIONAL  PHASES  SHOW  ANDEAN  ZONE  COMPLEXITY 

S.  Minaya,  R.  Cabre  S.  J.  and  R.  Ayala 
Observatorio  San  Calixco.  Cas.  12656.  La  ?az,  Bolivia 


Regional  phases  Li,  Lg,  Rg  and  5,  ?  residuals,  for  Sourh  America!i 
earthquakes  recorded  at  L?3  station,  show  a  regional  complexity  rlaiming  for  a 
new  division  in  subregions  of  different  seismic  and  structural  characteristics. 

Eastern  South  America,  excepting  the  quite  stable  Guyana  Shield,  has  a  few 
anomalies,  probably  originated  in  deep  fractures.  In  northern  part  or  South 
America  attenuation  is  higher  t.han  normal,  uneven,  wit.h  a  tendency  to  increase 
fio.-Ti  W  to  S.  The  Andes  in  ?eru-2cuador  are  the  most  complex  part  of  the  whole 
Scuti-.  .America,  especially  Lg  and  Rg  are  different  for  earthquakes  supposed 
si.milar  in  location  and  site.  Continuing  to  t.he  S,  we  find  double  pnase  ?,  but 
only  wnen  quakes  originate  i.n  suoducted  plate;  attenuation  of  Lg  still  appear:- 
nigh.  Scut.nerr.  Argenti.na-Ch.le  is  t.te  less  tom.plex  part  of  the  Andes,  though  1^., 
Lg  a.td  Rg  are  not  recorded  at  ail;  eventually  ?  phase  is  double  or  also  triple. 


U 


Fig. 


Ic.  Particle  motion  of  phase. 


.  PARTICLE  MOTION  OF  THREE  P  PHASES.  ZONE  8f  (Fig.  6.8). 


15 


i 


H,... 


”  Best  Available  Copy 


Fi>;ure  A.  Xvizcn  plate  and  depth  of  Moho 
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Li  OF  PERU  EARTHQUAKES  RECORDED  AT  LPB 

The  analysis  of  Li  waves  originated  in  Peru,  recorded  at  LPB 
(La  Paz-Bolivia)  station  continues  a  previous  study:  La  onda  Lg  a 
v.raves  de  los  Andes  y  registradas  en  la  estacion  de  LPB  (Alcocer, 
1989).  Generally  the  beginning  of  Li  waves  is  enough  clear. 

350  Events  occurred  in  Peru  from  1974  to  1986  could  be  useful 
for  that  analysis;  among  them  95  were  selected,  excluding  the 
smaller  ones,  several  superposed  with  other  signals,  etc.  On  the 
other  hand  two  recent  events  clearly  recorded  in  local  digital 
stations  were  added  for  comparaison  reasons  (Fig.  7) . 


Fig.  7.  Epicenters  of  earthquakes  studied  here 
(1974-1986) . 
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The  guided  waves  Li  were  found  by  M.  Bath  (1956),  meanwhile 
studying  in  detail  the  following  phases  Lg  and  Rg  (  Press  and 
Ewing,  1952) ;  they  propagate  along  a  low  velocity  layer,  supposed 
to  extend  between  Conrad  and  Mohorovicic  discontinuities;  that 
supposition  originates  in  the  analogy  of  those  waves  with  3b  waves 
observed  for  near  earthquakes . If  the  path  focus-station  is  ocenaic 
along  about  100  km,  Li  is  not  recorded,  but  when  the  oceanic  path 
is  no  more  than  60  km,  Li  is  recorded  the  same  as  when  the  path  is 
completely  continental.  The  non-existence  of  a  granitic  layer, 
together  with  the  reduction  of  basaltic  layer  to  a  few  km,  mean 
that  the  continental  crust  is  responsible  of  Li  generation  and 
transmission;  so,  a  low  velocity  layer  is  thought  to  occupy  the 
intermediate  space  (so  it  is  represented  by  index  i) . 

When  the  path  of  Li  waves  is  parallel  or  subparallel  to  .\ndes 
Cordillera,  Li  are  smaller  than  Lg,  being  the  ratio  Li/Lg 
0.65±0.20;  both  trains  of  waves  have  a  similar  predominant  period 
(that  allows  a  direct  comparaison  of  amplitudes  in  the  record 
itself) .  If  path  is  perpendicular  to  the  Cordillera,  Li  at  least 
equals  Lg  amplitude,  but  generally  is  larger,  being  the  mean  ratio 
1.05±0.21. 

Li  are  the  beginning  for  several  trains  of  surface,  short 
period,  low  velocity  waves,  appearing  as  a  continuation;  so,  some 
times  the  whole  set  of  different  waves  is  called  Lg. 
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Particle  motion  of  Li  waves  appears  complex,  mainly  for 
surface  earthquakes,  what  suggests  that  such  waves  are  a 
composition  of  several  arrivals. 

The  largest  part  of  the  energy  belongs  to  the  horizontal 
proiecrion  SH  (transverse  to  the  ray  path) ;  after  that  the  largest 
amplitude  corresponds  to  the  component  SVH  (parallel  to  the  ray 
path);  the  vertical  amplitude  SW  is  small,  practically  absent  for 
intermediate  depth  earthquakes. 

Amplitude  spectrum  was  obtained  by  fast  Fourier  transform; 
the  maximum  amplitude  corresponds  to  frequencies  o.75  to  1.40  Hz. 

WAVES . 

The  ground  apparent  velocity  of  Li  was  calculated  by  dividing 
the  hypocentral  distance  (D)  by  the  travel  time  (t) : 

V  =  D/t 

Two  different  ranges  of  velocity  may  be  distinguished 
corresponding  to  different  ranges  of  depth  (see  Fig.  8) . 


Depth  (km) 

Velocity  (km/s) 

hcl20 

h>120 

3 . 84±0. 03 

3.91+0.04 

7  7 


Fig.  8.  Histogram  showing  the  correspondence 
between  frequency  of  apparent  velocity  values 
considering  two  depth  ranges. 


23 


25 


PART  II 

ATTENUATION  OF  THE  SHORT  PERIOD  HAVES  P,  Li,  Lg  AND  Rg 
TRH0U6H  PERU  -  BOLIVIA,  RECORDED  AT  LPB 

R.  Rodolfo  S.  Ayala,  Ramon  P.  Cabre  S.J.  and 
Estela  R.  Minaya 

Observatorio  San  Calixto,  La  Paz,  Bolivia 


ABSTRACT 

This  preliminary  study  analyses  the  characteristics  of  phases  P, 
Li,  Lg  and  Rg  of  earthquakes  originated  in  Peru,  recorded  at  La 
Paz  -  Bolivia  (LPB).  See  the  mean  apparent  velocity  (km/s)  for  those 
phases,  according  the  depth: 


P 

Li 

Lg 

Rg 

h:S70  km 

7.7310.01 

3.8410.0001 

3.5410.0003 

3.1710.001 

70<hS217kra 

7.8610.05 

3.8110.01 

3.5310.003 

3.1510.00! 

The  anelastic  attenuation  has  values  quite  different  in  two 
separate  areas:  In  the  Cordillera  region  oscillates  between  0.12 
and  0.17  in  Brazilian  shield  and  Subandcan  oscillates 
between  0.12  and  0.17  Being  quality  factor  Q  inversely 
related  to  :In  the  Cordillera  region  Q  changes  between  192  and 
b30;  in  Brazilian  Shield  and  Subandcan  between  450  and  794.  The 
complexity  of  Cordilleran  structure  may  explain  a  stronger 
attenuation. 


Remark.-  The  Spanish  version  of  this  paper  is  in  press,  in  the 
Revista  Geofisica  (Pan  .hmerican  Institute  of  Geography  and  History). 
It  is  annexed  to  this  Report.  Figures  in  the  Report  are  taken  from 
the  Spanish  version. 
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INTRODUCTION 

Longitudinal  waves  (  P  waves)  have  particle  motion  in  the 
same  direction  as  the  waves  propagate.  Their  velocity  is  less  than 
3  km/s  in  the  crust  and  3  to  10  km/s  in  the  mantle.  They  are  well 
observed  in  the  short  period  vertical  component. 

Li  waves  (identified  by  Bath,  1956)  are  short  period  c.hannel 
waves  period  transmitted  along  the  intermediate  layer  of  the  crust. 
Their  group  velocity  is  3.79±0.07  km/s.  They  are  identified  as  S* 
wave  at  short  epicentral  distance,  t.heir  particle  motion  is 
complex,  with  large  horizontal  component. 

Lg  guided  waves  have  short  period,  large  amplitude  (discovered 
by  Ewing  and  Press,  1952)  ;  they  have  a  particle  motion 
perpendicular  to  the  ray  path,  polarized  on  horizontal  plane;  their 
group  velocity  is  3.54  km/s. 

Rg  were  recognised  (by  Bath,  1954)  as  guided  waves,  short 
period,  Rayleigh  wave  type,  with  propagation  along  the  granitic 
layer,  group  velocity  about  3.0±0.7  km/s  and  particle  motion 
characteristic  of  Rayleigh  waves. 

DATA 

The  events  studied  occurred  in  Peru,  2'’  to  18°  latitude  sourh, 
63°  to  81°  longitude  west  (Fig.  1) ,  magnitude  mb  4  to 
5.6,  depth  divided  in  two  ranges  h<70  km  and  70<h<2l7  km,  during 
1973  to  1986.  The  data  was  taken  from  ISC,  for  earthquakes  recorded 
at  LPB  station. 
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Pig.  1.  Epicenter  distribution  Map  according  to  depth. 
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METHOD 

The  maximum  amplitude  and  corresponding  period  of  P  waves  were 
measured  on  the  short  period  vertical  seismograms;  Li,  Lg  and  Rg 
maximum  amplitude  and  period  were  measured  on  short  period 
horizontal  seismograms. 

The  region  was  divided  in  six  diferent  zones  (Fig.  2) 
according  to  the  back  azimuth  station-epicenter.  In  relation  to 
the  depth  each  zone  has  two  levels  (Table  1) .  The  first  level  for 
shallow  depth  and  the  second  one  for  intermediate  depth. 


Amplitude  data  were  normalized  at  mb=5.0  in  order  to  avoid  the 


source  effect. 


fl.yv.  73V 


TABLE 

1 

1 

1 

j 

1 

ZONE 

AZIMUTH 

(sta-epi)  j 

i 

1 

340® 

1 

1 

-  32C°  ■ 

1 

2 

320° 

-  310°  : 

3 

310° 

-  280°  j 

280° 

-  27C°  1 

5 

270° 

-  250° 

6 

250° 

i 

-  160° 

1 

Fig.  2.  Epicenter  distribution  map, 

showing  the  division  by  zones. 
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Nuttli's  equation  of  magnitude  (1973)  was  used  to  compute  the 
C,  3  coefficents  for  each  phase  and  zone.  Coefficent  of  attenuation 
was  computed  from  attenuation  curve  (Fig.  3)  of  Nuttli  (1973), 
whose  basic  equation  is: 

A-  K*(l/  (0“)  )  '''^*(1/  (SIN(D“)  )  )  (6'^*°^ 

Where  A  is  amplitude  (micrometers),  K  constant,  D”  epicentral 
distance  and  T  anelastic  attenuation  coefficent  (’/degree). 


Pig.  3.  Attenuation  theorical  curves  of  the  amplitude 
in  time  domain,  where  <  is  expressed  in  reverse  of 
degree  (Nuttli,  1973)  . 
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Units  of  y  were  changed  to  ’/km  in  order  to  get  the  quality 
factor  Q  using  Nuttli  equation  (1973)  : 

Q  =7r/(T*i*V) 

Where  T  is  period  in  seconds,  ^  anelastic  coefficent  of 
attenuation  (’/km)  and  V  group  velocity  (km/s) 

Herrmann  (1980)  method  was  applied  to  data  in  order  to  verify 
the  Q  values  obtained  by  Nuttli 's  method  (1973).  Herrmann's  method 
used  the  coda  of  the  seismic  event,  considering  that  the  main 
frequency  fp  (Hz)  is  inversely  proportional  to  Q  values. 

The  master  curves  were  used  to  obtain  the  Q  values  (Fig.  4) 
for  short  period  standard  seismograph  (WWSSN) .  Travel  time  t 
divided  by  Q  is  equal  to  t*  (Fig.  4) ;  fp  (Hz)  was  computed  taking 
a  time  band-window  about  10  seconds  and  number  of  zeros  t.hat  cross 
the  mean  trace  of  signal,  after  that  they  are  divided  twice  by  the 
time  band-window;  fp  (Hz)  is  plotted  versus  t  (s)  and  cheking  wit.h 
the  master  curve  until  t*  becomes  equal  to  one,  the.n  the  Q  value 
will  be  numerically  equal  to  t. 


29 


31 
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tt  (s) 


Pig.  4.  Master  curve  of  main  frequency  fp  (Hz)  versus  t* 
for  standard  short  period  seismographs  (WWSSN)  (Herrmann, 
1980)  . 


Group  velocity  was  obtained  by  plotting  travel  time  versus 
epicentral  distance  fitting  curve  and  computing  it  for  each  phase 
(Figs.  5a  and  5b,  according  to  depth). 
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TECTONIC  AND  GEOLOGICAL  ENVIRONMENT 

The  region  studied  corresponds  to  a  convergency  zone  between 
two  tectonic  plates  where  the  Nazca  plate  is  subducted  under  the 
South  America  plate;  the  Nazca  plate  moves  a  6  cm/ year  and  the 
South  America  plate  3  cm/year  (Minster  and  Jordan,  1978)  . 

The  Nazca  plate  is  divided  in  five  segments  in  this  region, 
the  first  between  2“  and  6“  south  latitude  ,  dipping  12 °E,  100  km 
thick,  seismically  active  until  250  km;  the  second  segment  from  7 
to  12®  south  latitude  ,  dipping  20®E,  110  Km  thick  and  seismically 
active  until  580  km;  the  third  segment,  12®  to  15°  south  latitude 
,  dipping  16®E,  100  km  thick  and  seismically  active  until  about  290 
km;  the  fourth  is  located  from  15®  to  17°  south  latitude,  dipping 
21 '’E,  120  km  thick  and  seismically  active  until  about  300  km;  the 
fifth,  from  17°  to  24®  south  latitude,  dipping  29°E,  120  km  thick 
seismically  active  until  about  580  to  600  km  (Ayala,  1991). 

Both  in  the  Northern  and  Sourhern  parts  there  are  active 
volcanoes,  being  weakness  zones. 

Continental  plate  is  divided  (at  surface)  in  five  geological 
provinces:  Andean  Cordillera,  Altiplano,  Coastal  Cordillera, 

Coastal  Plateau,  Upper  Basin  of  Amazonas,  Chaco-Benian  plains, 
Brazil  and  Guayana  Shields. 
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Andean  Cordillera  consists  of  Western  and  Eastern  Cordillera, 
the  last  one  composed  of  Paleozoic  rocks  and  granitic  bodies,  the 
spurs,  called  Subandean,  are  composed  mainly  by  sedimentary 
Paleozoic  rocks  ,  strongly  folded  and  fractured.  Western  Cordillera 
consists  mainly  of  volcanic  rocks. 

Coastal  Cordillera  is  made  of  a  series  of  batholithic  rocks 
parallel  to  the  coast  line.  It  is  a  local  range. 

Coastal  Plateau  separates  the  Coastal  Cordillera  from  the 
Western  Cordillera. 

High  Plateau  corresponds  to  a  tectonic  trench  lifted,  of 
mesozoic  age.  It  is  located  between  Western  and  Eastern 
Cordilleras,  filled  with  sedimentary  rocks,  paleozoic,  cretaceous, 
tertiary  and  quaternary  cover,  around  12000  meters  thick  (Schlater 
and  Mederlof ,  1966)  . 

Brazilian  Shield  consists  of  precambrian  rocks.  They  are  the 
basement  of  South  American  plate.  The  part  outcroping  in  Bolivia 
is  labeled  Guapore  Shield;  the  northern  prolongation  is  labeled 
Guayana  Shield,  These  are  separated  by  sedimentary  paleozoic  rocks, 
of  marine  origin  and  meso-cenozoic  rocks  of  Amazonas  basin  700 
meters  thick  (Jeffreys  and  Aires,  1988). 
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Pig.  6.  Tectonic  geological  map  (adopted  from  the  Mapa  Noetectoncio 
Preliminar  para  Sudamerica,  1985) , 
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The  Eastern  part  is  called  Upper  Amazonas  Basin,  with 
Paleocene  rocks.  It  covers  the  contact  between  Shield  and 
Cordillera.  It  is  about  4000  meters  thick  (Morrison,  1962). 

Chaco-Benian  plains  are  wide,  with  Cambrian  and  paleozoic 
rocks,  2300  meters  thick  (Morrison,  1962) ,  while  along  the  Beni 
river  the  precambrian  basement  is  covered  with  terciary  rocks. 


SEISMIC  HAVES  ORIGINATED  IN  THE  PERU 

The  Fig.  7  shows  the  short  period  vertical  records  of 
Peruvian  earthquakes,  recorded  at  LPB. 


Fig.  7.  Vertical  short  period  seismograms  recorded  at  LPB. 
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(4) 
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km 
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P  Phase 


In  general  the  Peruvian  P  waves  have  emergent  beginnings: 
periods  between  0.6  and  1.2  s;  the  average  apparent  velocity  is 
7.73±0.01  km/s,  normalized  amplitudes  from  0.01  to  0.3  micrometers 
for  shallow  focus;  for  intermediate  focus  the  velocity  is  7.86±0.05 
km/s  and  their  normalized  amplitudes  from  0.02  to  i.l  micrometers. 
The  particle  motion  parallels  the  ray  propagation  (in  the  three 
planes:  horizontal  plane  formed  by  N-S  and  E-W  components;  vertical 
plane  parallel  to  the  azimuth,  formed  by  vertical  component  Z  with 
the  radial  proyection  UR,  and  the  vertical  plane  perpendicular  to 
the  azimuth,  formed  by  the  Z  component  and  the  horizontal 
proyection  UH  (Fig.  8a  y  8b) .  These  waves  travel  across  the  upper 
mantle,  for  the  considered  distances. 

Li  Phase 

They  are  guided  waves  of  short  period,  generally  with  emergent 
beginnings,  transmitted  by  the  lowest  part  of  the  continental 
crust;  periods  between  0.7  to  1.2  s;  average  velocity  is 
3.84t0.0001  km/s,  their  maximum  amplitudes  from  0.04  to  0.3 
micrometers  for  shallow  focus;  apparent  velocity  for  intermediate 
depth  is  3.81C0.01  km/s  and  their  maximum  amplitudes  from  0.1  to 
3,4  micrometers.  Particle  motion  shows  predominant  movement  almost 
perpendicular  to  the  ray  path  with  some  horizontal  polarization 
(Fig.  9a), record  is  clearer  for  intermediate  focus  (Fig.  9b);  this 
phase  may  be  a  superposition  of  higher  modes  of  Love  waves. 


6a  .  Movi  in  lento  de  parti'cuia  Tase  P, 


sismo  5upi?r  f 


Movimiento  .ie  r-artf.:ul-i 


aismu  inter medio 
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Lg  Phase 

They  are  guided  waves  of  short  period,  generally  with  emergent 
beginnings  ;  transmitted  by  the  higher  part  of  the  continental 
crust,  with  periods  between  0.6  to  1.5  s;  average  velocity  is 
3.54±0.0003  km/s,  their  maximum  amplitudes  from  0.02  to  0.9 
micrometers;  for  intermediate  focus  is  3.53±0.003  km/s  and  their 
amplitudes  from  0.09  to  2  micrometers.  Their  particle  motions  are 
perpendiculars  to  the  ray  path,  with  a  strong  horizontal 
polarization  (Fig.  10a) ,  major  in  the  intermediates  (Fig.  lOb) . 


Rg  Phase 

They  are  guided  waves  of  short  period,  generally  with 
emergent  beinnings;  interfered  by  the  previous  phases  coda;  it  is 
transmitted  by  the  higher  part  of  the  continental  crust; 
predominant  periods  between  0.6  to  1.4  s;  average  velocity  is 
3.17r0.001  km/s,  their  maximum  amplitudes  from  0.03  to  0.2 
micrometers  for  shallow  focus;  for  intermediate  focus  velocity  is 
3.15±0.005  km/s  and  their  maximum  amplitudes  between  0.04  to  1.1 
micrometers.  Particle  motion  is  parallel  to  the  ray  path  (Figs, 
lla  and  11b) ,  being  a  composition  of  much  higher  modes  of  Rayleigh 
phase. 
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ANELASTIC  ATTENUATION  COEFFICIENT 

In  order  to  estimate  the  anelastic  attenuation  coeficients 
(i'),  first  the  amplitudes  read  for  each  zone  and  for  differeni: 
depth  levels  were  normalizated,  to  mb=5.0.  Real  values  were 
superposed  to  the  theorical  curves  (Fig.  3)  obtaining  ^  the 
values  (in  1 /degrees)  for  each  considered  phase. 

The  Figs.  12a  to  18a  show  the  i  values  for  each  phase,  zone 
and  depth  level  analysed.  The  values  of  period  (T)  and  average  ^ 
and  their  respectives  standard  deviations  are  detailed  in  the 
Table  2. 

QUALITY  FACTOR 

With  the  obtained  values  the  formuled  (II)  was  applyed  to 
calculate  the  apparent  Q  values  for  each  zone,  considering 
different  focal  dept.hs. 

Plotting  fp  versus  t,  was  appropopiated  through  the  metod  of 
seismic  wave  coda  form  (Herrmann,  1980)  ,  that  may  confirm  the  Q 
values  calculated  by  the  Nuttli's  method  (Figs.  19a  to  26b). 

The  computed  Q  values  were  plotted  on  a  map  to  obtain  the 
regions  of  similar  Q  value  (Figs.  27a  to  30b)  . 
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Group  velocity  values  (V)  and  average  Q  for  each  zone  and  by 
depth  levels  and  their  standard  desviations  are  expressed  in  the 
Table  2. 


TABLE  2 


1  PHASE 

1  ZONE 

1  T 

(s) 

i  ^ 

1 1  /dearee » 

;  V 

( km/s ) 

:  <2n 

Qh 

1 

0.35±0.01 

7.7110.03 

2 

0.30+0.01 

0.1210.004 

7.8810.02 

45811 

450 

P 

■k 

3 

0.90±0.06 

7.8010.05 

4 

0.80+0.04 

0.1610.00007 

7.77+0.05 

357110 

350 

5 

0. 90±0.06 

0.30  0.005 

7.5910.02 

181  12 

200 

1 

0.35±0.02 

8.0110.03 

2 

0.35+0.04 

0.1210.0005 

7.9610.002 

447126 

450 

P 

n  «r 

3 

0.30±0.02 

7 . 9010.04 

5 

0.3010.007 

0.3010.0008 

7.6510.006 

192111 

200 

6 

0.9510.02 

6.9810.25 

1 

1.010.007 

3.3410.003 

2 

1.0 

0.1210.00005 

3.34+0.003 

75311 

7  50 

Li 

It 

3 

1.0 

3.3410.005 

4 

0.95  0.01 

3.34  0.007 

5 

0.9010.02 

0.1610.01 

3.3410.006 

60119 

600 

1 

0.9010.01 

3.3310.004 

2 

1.0 

0.1210.00008 

3.8110.004 

79415 

730 

Li 

k  * 

3 

1.0 

3.8110.001 

4 

1.0 

3.84 

5 

0.910.06 

0.1610.0002 

3.77*0.005 

630123 

600 

6 

1.0 

3.73 

1 

1.010.03 

3.5410.002 

2 

1.210.24 

0.1510.00004 

3.5410.004 

61418 

600 

Lg 

♦ 

3 

1.010.06 

3.5410.003 

4 

1.010.06 

3.54+0.001 

5 

1.010.03 

0.2010.006 

3.5410.003 

48015 

480 

1 

1.0 

3.54*0.001 

2 

1.010.03 

0.1510.00003 

3.54+0.004 

63015 

650 

Lg 

*  * 

3 

0.9010.04 

3.5410.002 

5 

1 . 0 

0.2010.0004 

3.5210.003 

51216 

500 

1.010.07 

0.2310.06 

3.4810.07 

420110 

400 

1 

1.010.03 

3.1810.004 

2 

1.010.01 

0.1610.003 

3.1810.004 

66414 

650 

Rg 

♦ 

3 

1.110.03 

3.18*0.006 

4 

1.0 

3.1810.003 

5 

1 . 0510.01 

0.2010.002 

3.17+0.004 

529119 

500 

6 

1.1010.07 

3.17+0.002 

1 

1.010.06 

3.1710.002 

o 

A. 

1.010.006 

0.1710.002 

3.1710.001 

652123 

650 

Rg 

*  ★ 

3 

1.010.10 

3.1610.003 

4 

1.010.06 

3.  mo.  01 

5 

1.110.04 

0.2010.0008 

3.1110.004 

53416 

500 

6 

1.010.03 

3.1110.009 

•  shallow  focus  **  intermediate  focus 
N  Nuttli's  method  (1973)  H  Herrmann's  method  (1980) 
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INTERPRETATION 


Comparing  the  group  velocity  calculated  for  the  P  waves  with 
the  regional  velocity  model  (Fig.  32)  (Ayala,  1991) ,  we  notice  that 
these  waves  are  transmitted  through  the  higher  mantle,  including 
its  low  velocity  layer.  For  shallow  focus  three  regions  are 
distinguished  in  relation  to  Q,  the  first  one  with  lowest 
attenuation  corresponds  to  Brazilian  Shield  and  Subandean  zone  with 
the  highest  values  of  Q;  the  second  zone  corresponds  to  the 
Cordillera,  with  a  larger  attenuation,  caused  without  doubt  by  the 
structure  complexity,  and  the  third  sone  for  the  Altiplano  and 
active  part  of  Western  Cordillera.  For  intermediate  focus  two 
regions  are  distinguised  with  an  increased  attenuation  towards  the 
West;  the  first  for  the  Brazilian  Shield,  Subandean  and  Eastern 
Cordillera  and  the  second  for  the  Altiplano  and  Western  Cordillera. 

The  Q  calculated  values  of  Table  2  correspond  to  an  average 
of  the  whole  trajectory;  analysing  these  values  as  function  of 
the  epicentral  distance,  a  tend  to  a  little  increase  of  Q  with  the 
distance  appears,  that  will  be  studied  larer.  For  more  distant 
events  there  is  a  deeper  penetration  of  the  P  waves  corresponding 
zo  increased  density  and  homogeneity  of  mantle  material. 


50 


53 


The  Fig.  3  3  shows  the  depth  variation  of  the  Mohorovicic 
discontinuity  between  Peru  and  Bolivia  (modified  from  James,  1971) , 
the  maximum  thickness  of  the  crust  below  the  central  Andes 
is  70  km  (Fernandez,  1968;  James,  1971  and  Molina,  1977),  it 

diminishes  eastwards  until  40  km,  to  become  constant  with  a 
thickness  of  39  km  for  the  Brazilian  Shield  region  (Oblitas,  1972) ; 
westwards  the  thickness  of  oceanic  crust  in  the  Nazca  plate  is  5 
kilometers  (Guzman,  1972) . 

Two  cases  of  guided  waves;  shallow  focus,  two  regions 
according  to  Q  values,  may  be  distinguished  concerning  attenuation 
is  small  in  the  shield  region,  it  increases  westwards  beneath  the 
Cordillera. 

For  intermediate  foci  there  is  a  better  resolution,  possibly 
because  the  attenuating  effect  of  the  sedimentary  package  of  the 
continental  crust  is  lesser;  studyng  the  quality  factor  values,  we 
can  discriminate  two  regions  and  also  three  for  Lg;  noticing  a 
larger  Q  value  (minor  attenuation)  for  the  Brazilian  Shield  and 
Shield-Cordillera  border  (Subandean) ;  Q  diminishes  in  the  Eastern 
Cordillera  (larger  attenuation),  and  much  more  in  the  Altiplano 
and  Western  Cordillera;  these  vertical  and  lateral  changes  in  the 
structure  appear  trough  Q  changes;  sea  it  in  a  preliminary  model 
of  Q  for  the  continental  crust  between  Peru  and  Bolivia  (Fig.  34) 
to  show  the  complex  structure  below  the  Central  Andes. 
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Fig.  34 


Preliminar  model  of  Q  for  the  continental  crust 
between  Peru  and  Bolivia  (according  to  James,  1976) 


The  attenuation  study  of  guided  waves  of  short  period  and 
intermediate  focus  allows  to  discriminate  better  different 
structures . 
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ANNEX  1 


The  previous  preliminary  study  was  proposed  for 
publication  in  Spanish  in  the  Revista  Geofisica  (Pan 
American  Institute  of  Geography  and  History)  and  is  in 
press  when  this  report  is  written.  It  is  annexed  to  the 
report . 
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ATENUACION  DE  LAS  FASES  DE  CORTQ  PERIODO  P,  Li,  Lg  y  Rg 
A  TRAVES  DE  PERU  -  BOLIVIA,  REGISTRADAS  EN  LPB 


Rene  Roaol to  Ayala  S. 
Ramon  Cabre  Roige  S.J. 
Esteia  ilinaya  Ramos  * 


ABSTRACT 

This  preliminary  study  analyzes  the  characteristics  of  phases  P, 
L;,  Lg,  and  Rg  of  earthquakes  originated  in  Peru,  recorded  at  La  Paz- 
Solivia  (LPB). 

P  L:  Lg  Rg 

<m  '.rCta.Jl  3.^4+0.3001  3.J4i0.3003  3.l7t0.30l 

2173h>70km  7.36i0.3S  2.3110.01  3.5310.002  3.1510.005 

The  anelastic  attenuation  ^  has  values  quite  different  in  two 
separate  areas:  In  Cordillera  region  s' osc  1 1 1  ates  between  0.16  and  3.30 
1/degree;  in  Brazilian  Shield  ana  Subandean  ^oscillates  between  3.12 
and  3.17  degree. 

3eing  quality  •' actor  3  inversely  '■elated  to  JT  :  In  the  Cordillera 
region  3  manges  between  '92  and  530:  in  Brazilian  Shield  ana  Subanaean 
oetween  450  and  7?4.  The  comolexity  qt  Cordilleran  structure  may 
axoiain  a  stronger  attenuation. 


3bser  vac  or  .  o  Ian  Calixto.  La  ^’az.  Bolivia 
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Rssijm. 


El  presente  estudio  prel iminar  analiza  las  car acter ist icas  de  las  'ases 
P,  Li,  Lg  y  Pg,  de  sismos  aei  Peru,  regtstradas  an  La  Paz  -  dolivia 
(LP9). 

La  veiocidad  promedio  (km/s;  nara  las  fases  csnsideradas  es: 

P  Li  Lg  Pg 

h<70  km  7.73ifl.ai  .  3.34+0.0001  3.54+0.0003  3.1710.301 

217>h>70km  7.36+0.03  3.3110.01  3.33+0.003  3.15+0.303 

Se  tisnen  an  general  dos  valores  del  coefic rente  de  atenuaci6n  anelis- 
ticaf,  para  la  regidn  de  la  Cordillera  S'»0.12  a  0.17  1/grados  y  para 
el  Escudo  Brasileno  y  Suoandino  ‘<'=0.16  a  3.30  1/graaos;  en  reiacron 
inverse  con  el  factor  de  calidad  '3;  sienoo  J»192  a  630  y  0*450  i  794 
respect ivamente;  corresoondiendo  la  mayor  atenuaci6n  de  las  ondas  a 
la  region  de  la  cordillera.  deOido  a  la  comolejidad  en  la  estructura. 


[NTRODUCCION 

Las  ondas  orimarias  P.  u  ondas  longitudinales.  tienen  un  movimiento  de 
partlcula  paralelo  a  la  direccidn  de  prooagacidn  del  rayo.  ^-'esentan 
velocidades  menores  a  3  km/s  en  la  corteza  y  de  3  a  10  km/s  en  el  manto.  Be 
registran  mejor  en  los  sismogramas  vert :c ales  de  corto  periodo. 

Las  ondas  Li  ( ident i f icadas  por  3ath,  1956)  cprresponoen  a  onaas  guiaaas 
de  corto  periodo.  que  se  transmicen  por  La  caoa  intermedia  at  la  corteza  Cde 
ahl  el  subindice  i),  con  velocidades  de  grupo  ae  3.79t0.37  km/s.  /  eouivalen 
a  la  fase  3*  observada  para  sismos  cercanos;  presentan  un  movimiento  oe 
particula  comolejo.  :on  una  mayor  comoonente  norizontal. 

Las  ondas  Lg  son  ondas  guiaUas  3e  corto  periodo  /  gran  amolitua.  pue  se 
suoone  se  orooagan  oor  la  oenominada  :aoa  gr-initica  de  la  corteza.  te  ani  el 
suomdice  g;  las  descuor ;  eron  Eving  /  ^’’ess  en  1952;  oresentan  un  movimiento 
de  oarticula  oeroend icui ar  a  la  direction  pe  ornoagacion  jei  <'avc.  oo-arizacas 
?n  ?i  piano  horicontal,  :on  /eiocidaa  oe  gruoo  de  1.-4  -.m/s.  ze  'egisc-'an 
nejor  en  los  sismogramas  norizontales  Pe  corto  per-ooo. 

Las  ondas  Pg  freconocidas  por  3ach.  I9;i4'i  jon  onoas  ;ui«aas  ;e  corto 
per iodo.  de  tipo  Payleiqn.  due  se  transmiten  oor  la  denomtnaaa  laca  ;raniCica 
de  la  corteza;  oresentan  un  movimiento  ae  jarticula  paralelo  a  la  direct  i<)n  de 
prooagacion  del  "ayo.  eliatico  /  retngraoo  en  el  oiano  /erticai;  una 
/elocioad  de  gruoo  de  3. 3+0. 7  vm/s.  Se  '•egistran  nejor  en  la  romoonente 
vertical  de  corto  oer-.odo. 
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La  regi6n  de  estudio  estA  comorendida  entre  los  2*  a  IB®  latitud  sur  y 
68®  a  81®  longitud  oeste;  se  consideraron  los  sismos  provenientes  del  Peru 
(Fig.  1),  con  magnitudes  mb  (magnitud  de  las  ondas  de  cuerpo)  desde  4.0  a  5.6, 
de  los  aflos  1977  a  1906,  koinados  de  los  Cak^logos  de  Sismos  del  Internat  lonal 
Seismological  Center  (I.S.C.),  registrados  en  las  componentes  de  corto  periodo 
de  la  eskaci6n  sismol6gica  de  La  Pat,  Bolivia  (LPB)  de  la  WWSSN  (World  Wide 
Seismographic  System  Network). 

Los  sismos  utili2ados  varian  desde  profundidades  superf ic lales  (h£70  km) 
hasta  interroedios  con  profundidades  (70<h<217  km). 


PROCEDIhlENIO  Y  METODO 


Para  medir  las  amplitudes  se  consideraron  los  siguientes  aspectos:  para 
las  ondas  P  se  midieron  las  m<liximas  amplitudes  y  los  periodos  correspond ientes 
de  los  sismogramas  verticales  de  corio  periodo  y  para  las  ondas  guiadas  Lg,  Li 
y  Rg,  las  m^ximas  amplitudes  y  periodos  de  los  sismogramas  horizontales  de 
corto  periodo. 

La  region  de  estudio  se  dividi6  en  6  zonas  diferentes  (Fig.  2),  de 
acuerdo  al  azimut  (estac idn-epicentro) .  Considerando  la  divisi6n  por 
profundidad  para  cada  zona,  se  tienen  dos  niveles  por  cada  una  de  el  las  (Tabla 
l)j  el  primero  con  los  sismos  super f ic lal es  y  el  segundo  con  los  intermedios. 


TABLA  1 

ZONA 

AZI11UT  (est-epi) 

1 

340®  -  320® 

2 

320®  -  310® 

3 

310®  -  280® 

4 

280®  -  270® 

5 

270®  -  250® 

6 

250®  -  160® 

4 


Con  los  datos  de  amplitud  y  periodos  leidos,  be  procedi6  a  normal iz ar  la 
amplitud  para  cada  zona  a  una  magnitud  mb=5.0,  considerando  similares  niveles 

de  profundidad,  para  evltar  el  efecto  de  foco;  a  partir  de  la  f6rniula  de 

magnitudes  (Nuttli,  lS/3): 

«b=  C  Bdog  Q«)  +  log(A/T)  (1) 

Oonde: 

mb  ~  maqniiud  de  las  ondas  de  cuerpo 
C  -  ufflbral  de  detec  tab 1 1 idad 
B  -  coeficiente  de  aienuacidn 
0*  s  distancia  epicentral  en  grados 
A  -  amplitud  maxima  eri  micrdmetros 

T  =>  periodu  de  las  ondas  de  amplitud  maxima  en  segundos 

Primero  se  determinarori  los  coef ic lenies  C  y  B  para  cada  fase  y  zona; 

despu^s  se  despejd  de  la  formula  (I)  la  amplitud;  reemplazando  todos  los  datos 
conocidos  se  normal iz6  la  amplitud  a  una  deierminada  magnitud  (mb=5.ei)  y  rango 
lia  nrn  fund  idad. 

Seguidamente  se  utilizaron  las  curvas  de  atenuaci6n  (Fig.  3)  de  Nuttli 
(lS/3)  para  calcuiar  el  coeficiente  de  atenuacidn  en  base  a  la  ecuacibn: 

V*  -^*0* 

K*(1/(C1«))*(1/(SIN(0«)))»(b  ) 


Oonde: 

A  =  amplitud  de  la  onda  en  micrbmelros 
K  =  constante 

0*  =  distancia  epicentral  en  grados 

If  coeficiente  de  atenuacibn  anel^stica  en  (l/grados) 


:U5rAI;CIA  Kl  irFMTHAI.  <J| 

Fig.  3.  Curvas  tebricas  de  atenuacibn  de  la  amplitud  en  el  dominio  del 
tiempo,  donde  S  estA  expresado  en  el  inverso  de  grades  (Nuttli, 
1973). 
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Para  calcular  el  factor  cJe  calidad  Q  se  camdiaron  las  unidades  de  a 
1/km  y  se  empled  la  formula  (Nuttli,  19/3): 


Q  =-n/(T«<*V)  (II) 


Oonde  i 

Ti  =3.  1416 

T  =  periodo  en  seguitdos 

'i  =  coeficiente  de  atenuaci6n  aneiastica  en  l/kil6nieiros 
V  =  velocidad  de  grupo  en  kil6metros  por  seguiido 


Para  comprobar  los  valores  de  Q  obtenidos  con  el  n^todo  de  Nuttli,  se 
utiliz6  tarbien  el  m^tudo  de  la  forma  de  la  coda  (Herrmann,  19B0),  el  cual 
explica  que  la  frecuencia  predominante  (fp)  es  inversamente  proporcional  a  □, 
para  lo  que  ut 11 izaremos  la  curva  maestra  fp  versus  (Fiq.  4)  para 
instrumentos  de  corto  periodo  de  la  red  standard  (WUBSN),  donde  t*  es  iqual 
t/Q,  slendo  t  el  tiempo  desde  el  urigen  en  segundos.  Para  este  propiisito  se 
calcula  la  fp  de  cada  fase  analizada  en  Hertz,  que  consiste  en  tomar  un  tiempo 
de  medida  10  segundos  o  mils  (ventana  de  tiempo)  y  contar  el  niimero  de  ceros 
que  cruzan  la  traza  media  del  sismoqrama  y  dividirlos  por  dos  veces  el  tiempo 
(ancho)  de  la  ventana;  se  plutean  fp  (Hz)  versus  t  (s)  y  son  sobreimpuestos 
encima  de  la  curva  maestra  (Fiq.  4)  hasla  que  l*-l,  entonces  el  valor  de  Q 
seri  num^r ic ament e  igual  a  t. 


t*  (s) 


Fiq.  4.  Curva  maestra  de  la  frecuencia  pr edominante  de  la  coda  fp  (Hz)  versus 
t»  (s)  para  instrumentos  de  corto  periodo  de  la  WWSSN  (Herrmann, 
1980). 


67 


6 


Para  determinar  la  velocidad  de  grupo  de  cada  fase,  se  plutearon  los 
tiempos  de  viaje  (tiempo  origen  iiit>nus  t  iiimpo  de  llegada  de  la  fase;  en 
segundos)  y  la  distancia  epicenlral  en  grades  (0*),  de  acuerdo  a  la 
profundidad  (Fig.  3a  y  3b>,  hactendu  un  ajuste  lineal;  para  despu^s  calcular 
la  velocidad  de  grupo  para  cada  tipo  de  onda. 


Tig.  3a.  Curvas  de  tiempo  viaje  versus  distancia  para  sismos  super f ic iales. 


Fig.  5b.  Curvas  de  tiempo  de  viaje  versus  distancia  para  sismos  intermedios. 
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Los  valores  del  factor  Ue  cal idad  Q  obtenidos  fueron  ploteados  en  un 
mapa  y  con  ayuda  de  la  tectdnica  y  geologia  de  la  regi6n  se  real iz6  la 
interpretac i6n. 


MARCU  fLCTONlCQ  -  GLOLQGICQ 


La  reqi6n  de  estudio  corresponde  a  una  zona  de  converqenc la  de  dos 
placas  fcect6nicas,  donde  la  piaca  de  Nazca  se  escA  subdue lendo  debajo  de  la 
piaca  Sudamer icana;  la  velocidad  de  avance  de  la  piaca  de  Nazca  es  de  Cj  cm/afio 
y  para  la  piaca  Sudamer  icaiia  de  3  cm/afio  (Hinster  y  Jordan,  1S78) . 


La  piaca  de  Nazca  en  la  region  de  estudio  presenta  cinco  segtnefitos;  el 
prltnero  entre  los  2*  a  6*  latitud  sur,  con  un  buzamiento  12*E,  espesor  de  100 
km  y  actividad  sismica  hasta  la  profundidad  de  2tj0  km;  el  sequndo  desde  los  7* 
a  12*  latitud  sur,  con  un  buzamiento  de  20*E,  espesor  de  110  km  y  profundidad 
de  500  km;  el  tercero  desde  los  12*  a  15*  latitud  sur,  con  un  buzamiento  de 
16*E,  espesor  aproximado  100  km  y  profundidad  de  290  km;  el  cuarto  desde  los 
15*  a  17*  latitud  sur,  buzamiento  de  21*E,  espesor  de  170  km  y  profundidad  de 
300  km;  el  quintu  desde  17*  a  24*  latitud  sur,  buzamiento  de  29*E,  espesor  do 
120  km  y  profundidad  de  500  a  600  km  (Ayala,  1991). 


En  la  parte 
de  debilidad. 


norte  y  sur  se  presentan  volcanos  activos,  mostrando  zonas 


En  la  piaca  continental  se  tienen  varias  provincias  qeoldqicas  (Fig.  6): 
Cordillera  de  los  Andes,  Altiplano;  Cordillera  de  la  Costa;  Planicie  Costera; 
Cuenca  Superior  del  Amazonas;  Llanuras  Chaco-Beruanas;  Escudo  Brasil eno  y  el 
CscuwO  ds  Quayana. 


La  cordillera  de  los  Andes  se  divide  en  cordillera  Oriental  y 
Occidental;  la  primera  const ituida  por  rocas  paleozoicas  y  muchos  cuerpos 
intrusivos  graniticos;  las  estr ibac lones  orientales  de  la  misma  se  denominan 
Subandino  y  son  rocas  sedimentar ias  paleozoicas,  fuertemente  pleqadas  y 
fracturadas;  la  cordillera  Occidental  esta  constituida  por  una  sene  de 
volcanes  y  mesetas  volcanicas. 


La  cordillera  de  la  Costa,  constituye  una  sene  de  elevaciones  paraleias 
a  la  costa  en  la  parte  sur  del  Peru,  formada  por  cuerpos  batoliticos;  sienUo 
una  cordillera  de  tipo  local. 

La  Planicie  Costera,  separa  la  cordillera  de  la  Costa  de  la  cordillera 
Occ idental . 

El  Altiplano  corresponde  a  una  fosa  tectbnica  levantada,  de  edad 
mesozoica  terciaria,  ubicada  entre  la  cordillera  Oriental  y  Occidental, 
rellenada  por  rocas  sedimentar  ias,  paleozoicas,  cret^cicas,  tercianas  y 
cubierta  cuaternaria,  con  un  espesor  de  12000  m  (Schlater  y  Nederlof,  1966). 
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El  Escudo  Brasileno  eslA  formado  por  rocas  de  edad  precdmbr ica,  que 
constituyen  el  basamento  del  continente  sudamer icano  (la  parte  que  aflora  en 
Bolivia  se  denomina  Crat6n  de  Guapor^);  la  prolonqaci^n  de  ^sie  en  la  parte 
norte  se  denomina  Escudo  de  Guayana;  separados  por  rocas  sedimentar las  de  edad 
paleo2oica  de  origen  marino  y  rocas  meso-cenozoicas  de  la  cuenca  del  Amazonas; 
^sta  tiene  un  espesor  de  7000  m  (Jeffreys  y  Aires,  1968);  la  parte  mas 
occidental  se  denomina  Cuenca  Superior  del  Amazonas,  con  rocas  paleo- 
cenozoicas,  que  cubren  el  contactu  entre  el  escudo  y  la  cordillera,  con  un 
espesor  aproximado  de  4000  m  (Morrison,  1962). 

Las  Llanuras  Chaco-Bentanas  sun  ampl ias  llanuras  const i tuidas  por  rocas 
paleozoicas  y  cambricas,  con  un  espesor  aproximado  de  23'10  m  (Morrison,  1962); 
mientras  que  a  lo  largo  del  rio  S^ni  el  basamento  precambrico  solo  esta 
cubierto  por  rocas  terciarias. 


UNDAS  SI SMI CAS  URIGINADAS  EN  EL  PERU 


La  Fig.  7  muestra  los  registros  verticales  de  corto  periodo  para  sismos 
del  Peru,  registrados  en  LPB. 


(>  ,  i.  1-1 

I  I  _  _ _ > _ i  1. 
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Fig.  7.  Sisfflogramas  verticales  de  corto  periodo  registrados  en  LPB. 
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Faae  P 


La&  ondas  P  dtfl  Peril  eri  general  tjenen  comienzo  emerqenie;  presentan 
periodos  entre  los  0.6  a  1.2  s;  su  velocidad  promedio  es  7.73i0.01  km/s,  sus 
amplitudes  normalizadas  variaii  entre  lus  0.01  a  0.3  micr4metros  (1 
micr6metro=0. 001  mm)  para  focus  super f icial es;  para  los  intermedios  es 
7.86i0.0t>  km/s  y  sus  amplitudes  varian  entre  0.02  a  1.1  micrdmetros.  Exhiben 
un  movimiento  de  particula  paralelo  a  la  direcci6n  de  propagaci6n  del  rayo  (en 
los  tres  pianos:  piano  horizontal  formado  por  las  componentes  N-S  y  E-W;  piano 
vertical  paralelo  al  azimut,  formado  pur  la  componente  vertical  Z  con  la 
proyecci6n  radial  UR  y  el  piano  vertical  perpendicular  al  azimut,  formado  por 
la  componente  2  y  la  pruyecci6n  horizontal  UH  (Figs.  8a  y  8b).  Son  ondas  gue 
viajan  por  la  parte  superior  del  manto,  para  las  distancias  gue  en  nuestros 
casos  interesan. 


Fase  Li 


Son  ondas  guiadas  de  corto  periodo,  con  comienzo  general merite  emerqente, 
gue  se  transmiten  por  la  parte  inferior  de  la  corteza  continental;  sus 
periodos  est^n  entre  los  0.7  a  1.2  s;  su  velocidad  promedio  es  de  3.8410.0001 
km/s,  sus  amplitudes  varian  entre  los  0.04  a  0.8  micr6metros  para  focos 
super f ic iales;  para  intermedios  es  3.8110.01  km/s  y  sus  amplitudes  entre  0.1  a 
3.4  micrbmetros.  Presentan  un  movimiento  de  particula  pr6ximo  a  la 
perpendicular  a  la  direccidn  de  propaqaci6n  del  rayo  con  alguna  polarizaci6n 
horizontal  (Fig.  9a),  mJis  claro  para  los  focos  intermedios  (Fig.  9b);  pudiendo 
ser  una  sobreposic 1 6n  de  modos  super  lores  de  la  onda  Love. 


Fase  Lg 

Son  ondas  guiadas  de  corto  periodo,  generalmente  Claras,  pero  con 
comienzos  no  impulsivos;  se  transmiten  por  la  parte  superior  de  la  corteza 
continental,  con  periodos  de  0.6  a  l.b  segundos;  su  velocidad  promedio  es 
3.J410.0003  km/s,  sus  amplitudes  varian  de  0.02  a  0.9  micrdmetros;  para  focos 
intermedios  es  3.b310.003  km/s  y  sus  amplitudes  entre  0.09  a  2  micr 6metrc»s. 
Sus  movimientos  de  particula  son  perpendiculares  a  la  direccidn  de  propagaci6n 
del  rayo,  con  una  fuerte  polarizaci6n  horizontal  (Fig.  10a),  mayor  en  los 
intermedios  (Fig.  I0b);  corresponded  a  modos  super  lores  de  la  fase  Love. 


F  ase  Rg 


Son  ondas  guiadas  de  corto  periodo,  con  comienzos  generalmente 
emergentes;  inter fer Idas  por  la  coda  de  las  fases  previas;  gue  se  propagan  por 
la  parte  superior  de  la  corteza  terrestre;  presenla  periodos  entre  los  0.6  a 
1.4  s;  su  velocidad  promedio  es  de  3.1710.001  km/s,  sus  amplitudes  varian 
entre  los  0.03  a  0.2  micrbmetros  para  focos  super f ic lal es;  para  focos 
intermedios  es  3.rji0.00'j  km/s  y  sus  amplitudes  entre  0.04  a  l.l  micrdmetros. 
Muestran  un  movimiento  de  particula  paralelo  a  la  direcci6n  de  propagacibn  del 
rayo  (Figs.  11a  y  lib);  siendo  un  composicidn  de  varios  modos  superiores  de  la 
fase  de  Rayleigh. 
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Fiy.  11a  •  Movimi^?nto  de  parti.;ula  I- asj  Pg,  sismo  sup^ir  f  ic  lal 


Fig.  11b  .  Movimiento  de  parti<;ula  F  ase  Pq,  sismo  intermedio 
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Para  determiriar  los  coeficientes  de  alE>nuaci6n  anel^Blica  (."!),  urimero 
se  normali^aron  las  amplitudes  leidas  para  cada  zona  y  por  niveles  de 
profundidad,  a  fflb=b.l9.  Ajustando  los  valores  a  las  curvas  te6ricas  (Fig.  3), 
obteniendo  los  valores  de  ^  (en  1/grados)  para  cada  fase  consider ada. 

Las  Figs.  12a  a  19b  muestran  los  valores  de  H  para  cada  fase.  zona  y 
niveles  de  profundidad  analizados;  los  valores  del  periodo  (T)  y  X  promedios  y 
sus  respect ivas  desviaciones  standard  esldn  detail  ados  en  la  Tabla  2. 


FACTOR  DE  CAL I DAD 


Con  los  valores  de  IT  ubterudos  se  aplic6  la  formula  (II)  para  calcular 
los  valores  de  Q  aparente  para  cada  zona,  considerando  diferentes  niveles  de 
profundidad  de  foco. 

Ploteando  fp  y  t,  se  procedi6  a  determinar  Q  por  el  m^todo  de  la  forma 
de  la  coda  de  las  ondas  sismicas  (Herrmann,  1980),  que  serviran  para 
corroborar  los  anter lores  valores  calculados  de  Q  (Figs.  20a  27b). 


Los  valores  de  Q  obtenidos  fueron  ploteados  sobre  un  mapa  para  obtener 
regiones  de  similar  valor  de  Q  (Fig.  28a  31b>. 


Los  valores  de  velocidad  de  grupo  (V)  y  0  promedios  para  cada  fas3  y  por 
niveles  de  profundidad  y  sus  desviacioties  standard  estan  expresados  en  la 
Table  2. 
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INTERPRETACIOW 


Conparando  la  velocidad  de  grupo  calculada  para  las  ondas  P  con  el 
modelo  de  velocidad  (Fig.  32)  para  la  reqidn  (Ayala,  1991),  notanos  que  $stas 
se  transfliiten  hasta  el  manto  superior,  incluyendo  su  capa  de  baja  velocidad. 
Para  locos  super ficiales  se  distinguen  tres  reqiones  en  relaci6n  a  Q,  la 
prifliera  con  menor  atenuaci6n  correspondiente  al  Escudo  Brasileno  y  Subandino, 
con  los  valores  mis  altos  de  Q;  la  segunda  2ona  corresponde  a  la  cordillera, 
con  una  mayor  atenuacibn,  debido  sin  duda  a  la  complejidad  en  la  estructura.  y 
la  tercera  para  el  Altiplano  y  parte  activa  de  la  cordillera  Occidental.  Para 
focos  intermedios,  se  presentan  dos  reqiones  con  un  aumento  de  la  atenuacidn 
en  direcci<in  oeste;  la  prtmera  para  el  Escudo  Brasileno,  Subandino  y  parte 
cordillera  Oriental  y  la  segunda  para  el  Altiplano  y  cordillera  Occidental. 


Los  valores  calculados  de  Q  de  la  Tabla  2  corresponden  a  un  proaedio  de 
toda  la  trayectoria;  anal iz^ndolos  en  funci6n  de  la  distancia  epicentral,  se 
presenta  una  tendencia  a  un  1  igero  auiaento  de  Q  con  la  distancia,  que  se 
estudiard  poster lormente.  Para  los  eventos  m4s  distantes  hay  una  mayor 
penetraci6n  de  las  ondas  P  y  un  consiguiente  aumento  en  la  densidad  y 
hoffiogeneidad  del  material  del  manto. 


La  Fig.  33  muestra  la  variacidn  de  la  profundidad  de  la  discont inuidad 
de  (lohorovicic  entre  PerCi  y  Bolivia  (mod.  de  James,  1971),  el  espesor  m^ximo 
de  la  corteza  debajo  de  los  Andes  centrales  es  de  70  km  (Fernindez,  1966, 
James,  1971  y  Molina,  1977),  el  cual  disminuye  en  direcci6n  este  hasta  40  km, 
para  mantenerse  casi  constante  con  un  espesor  de  39  km  para  la  regi6n  del 
Escudo  Brasileno  (Oblitas,  1972);  hacia  el  oeste  el  espesor  de  la  corteza 
oce^nica  en  la  Placa  de  Nazca  es  de  5  kildmetros  (GuzmJin,  1972). 


Para  las  ondas  guiadas  se  tienen  dos  casos;  para  los  focos  superf iciales 
se  presentan  dos  reqiones  de  valores  de  Q,  con  una  menor  atenuaci<)n  en  la 
regi6n  del  escudo,  la  cual  aumenta  en  direcci^n  oeste  en  la  cordillera. 


Para  los  focos  intermedios  se  tiene  mejor  resoluci()n,  posiblemente 
debido  a  que  el  efecto  atenuador  del  paquete  sedimentario  de  la  corteza 
continental  es  menor;  estudiando  los  valores  del  factor  de  calidad,  notamos 
que  se  pueden  discriminar  dos  reqiones  y  s41o  para  las  Lq  tres;  notando  un 
mayor  valor  de  0  (menor  atenuacidn)  para  la  reqi6n  del  Escudo  Brasileno  y 
tambi^n  para  el  limite  escudo-cordillera  (Subandino);  Q  disminuye  en  la 
cordillera  Oriental,  mayor  atenuaci^n,  y  aiin  aumenta  m^s  ^lla  en  el  Altiplano 
y  en  la  cordillera  Occidental;  estas  variaciones  verticales  y  laterales  en  la 
estructura  se  manifiestan  en  la  variacidn  de  Q,  expresado  en  un  modelo 
preliminar  de  Q  para  la  corteza  continental  entre  Peru  y  Bolivia  (Fig.  34), 
que  nos  muestra  la  complejidad  de  la  estructura  debajo  de  los  andes  centrales. 
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Fiq.  31  .  Modelo  prelimnar  de  Q  para  la  corteza  continental  entre  el  Peru  y 

Bolivia  (en  base  a  Jaaies,  1976). 


Por  todo  lo  anterior  riotaiaos  que  estudiando  la  atenuacidn  de  ondas 
gutadas  de  corto  periodo  y  focos  internedios  se  puede  discriminar  mejor 
diferentes  estructuras. 
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